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PACS. 71.70.Di - Landau levels. 

PACS. 73.20.-r - Electron states at surfaces and interfaces. 
PACS. 71.10.Ca - Electron gas. 



Abstract. - The statistical weight W as a, function of energy E for quasielectrons with mass 
m* subject to a fixed magnetic field B is W/A = (m* /nh 2 ) E + (eB/-K 2 ti) E^iM)^" 1 ' 
sin (2ni/E/hu!c), where A is the sample area, and u c = eB/m* the cyclotron frequency. Signif- 
icantly, there is no Landau's term proportional to B 2 in 2D. This leads to the conclusion that 
the 2D electron system is always paramagnetic, but shows a magnetic oscillation. 



Landau [1], Sondheimer and Wilson [2,3] discussed the de Haas- van Alphen (dHvA) oscil- 
lation [4] of a three-dimensional (3D) system of quasifree electrons. Nakamura [5] calculated 
the statistical weight W associated with the Landau states, and treated the dHvA oscillation. 
We extend Nakamura's theory to a 2D system. 

Let us consider a dilute system of electrons, each with effective mass m* , moving in a 
plane. Applying a magnetic field B perpendicular to the plane, each electron will be in a 
Landau state of energy 



where to c — eB/m* is the cyclotron frequency. The degeneracy of the Landau level (LL) is 



E = [Nx, + 1/2) hui c , N L = 0,1,--- 



(1) 



eBA 



A = sample area. 



(2) 



2ttH ' 



We introduce kinetic momenta 



n x = p x + eA. 



n 



y — Py + 



in terms of which the Hamiltonian TL for the electron is 




(3) 
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Fig. 1 - Quantization scheme for free electrons in 2D: without a magnetic field (dots), and in a 
magnetic field (circles). 



After simple calculations, we obtain 

dx dH x dy dH y = dx dp x dy dp y . 

We can then represent quantum states by quasi phase-space elements dxdH x dydH y . The 
Hamiltonian Ti. in eq. © does not depend on the position (x, y). Assuming large normalization 
lengths (Li, L 2 ), we can represent the Landau states by concentric shells of the phase space 
having statistical weig ht 27rII AIT • L 1 L 2 (2irh)- 2 = eBA/2nh, with A = L X L 2 and Huj c 
A (n 2 /2m*) = II AIT/m*. Hence, the LL degeneracy is given by eq. (J2J. Figure ^ represents 
a typical Landau state in the n^-IIy space. 

As the field B is raised the separation hoj c increases, and the quantum states are collected 
(or bunched) together. As a result of bunching, the density of states N{e) should change 
periodically since the Landau levels arc spaced equally. 

The electrons obey the Fermi-Dirac statistics. Considering a system of quasifree electrons, 
we define the Hclmholtz free energy T by 



T = Nfj, - 2k B T ■ In [l + e ( "- £i)/teT ] , (4) 



where /i is the chemical potential, and the factor 2 arises from spin degeneracy. The chemical 
potential /i is determined from the condition 

The magnetization M for the system can be found from 
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Equation JSJ is equivalent to the usual condition that the number of electrons N can be 
obtained in terms of the Fermi distribution function F: 

N = 2j2 F ( E i), F{E)=\e l3{E -^ + lY 1 . (7) 

i 

The LL Ei is characterized by the Landau oscillator quantum number {N-i) v 

Let us introduce the density of states dW/dE = AT(E) such that M{E)dE = number of 
states having energy between E and E + dE. We write eq. (@J in the form 

T = Nf,-2k B T- rdE^ln\l + e^/ k * T 
Jo dE I J 

poo 

= Nfi-2 / dE W(E)F(E). (8) 
Jo 

The statistical weight W is the total number of states having energies less than the Landau 
energy (N L + 1/2) hu c in eq. (|TJ). This W from fig. His 

L L °° 
W = j^^2ttUAU-2 ^ 6[£ - (iV L + 1/2) 7iw c ] , (9) 



where Q(x) is the Heaviside step function 

e(i) = 

We introduce a dimensionless variable e = 2itE/Tiu) c , and rewrite eq. © as 



if x < 

1 if x > 



VF(£) = C?iw c - 2^9[e- (2iV L + 1) tt] , (10) 

Af L =0 

with C = 2-Km* A(2nK)~ 2 . We assume a high Fermi-degeneracy such that /i ~ £f ^ ^w c . The 
sum in eq. (|10f> can be computed using Poisson's summation formula [6] 

OO -. oo />oo 

£ /(27m) = 5- E / d rf(T)e- imT , (11) 

n= — oo m— — oo 

where X)^L-oo /(2tt7i + i)j < t < 27r is by assumption uniformly convergent. We write the 
sum in eq. (fTUfl as 



2£e[e- (2n + ]>] = 6(e - tt) + 0), (12) 

oo 

^>(s;i) = E ©(e - tt - 27r|n + x|). 

Tl=— OO 

Note that <£(e; x) is periodic in x. After the Fourier expansion, we set x — to obtain eq. (|12fl. 
By taking the real part (Re) of eq. lfT21) and using eq. {TTJ, we obtain 



Re{ Eq. (HI} = - / dre(e-r) + -V(-l) l/ / drG(e - r) cos i/r, 

7f /q TT ^ Jo 



(13) 
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where we assumed e = 2irE/hij c 3> 1, and neglected 7r against s. The integral in the first 
term in eq. (|13[1 yields e. The integral in the second term can be evaluated by integration by 
parts, and using d<d/dy — S(y). We obtain 



1 

drQ(e — r) cos vt = — sin ve. 



Hence, 



1 2 °° f— lY 
Re{ Eq. O} = -e + - ^ - - sin ve. (14) 

77 y=l 



Using eqs. (fTTIj) and (fT4"j) . we obtain 

W(E) = W + W osc , (15) 

Wo = C7iw c ^ = (16) 
7r \ Tkcv / 



i/=i 



The oscillatory term Wosc contains an infinite sum with respect to v, but only the first term, 
v = 1, is important in practice as we see later. This term W osc can generate magnetic 
oscillations. There is no term proportional to B 2 generating the Landau diamagnetism. This 
is unexpected. We briefly discuss the difference between 2D and 3D systems. 
In 3D, the LL E is given by 

E = [ N h + 1} hw c ' ir 



2 J 2m* 

The energy E is continuous in the bulk limit. The statistical weight W' is the total number 
of states having energies less than E. The allowed values of p z are distributed over the range 

1/2 

in which \p z \ does not exceed {2m* [E — (N^ + 1/2) huj c }} 1 . We obtain 

/* x3/2 oo 

W'(E) = C' [ -^—2 £ y/e-(2N L + l)v, (18) 

where C = V (2irm*) 3 ^ 2 /(2?:h) 3 , e = 2 / nE/hoj Cl and V = L1L2L3 is the sample volume. We 
proceed similarly, and obtain 

W'(E) = Wq + W{j + W' oac , (19) 



W' Q = C'-^E 3 ' 2 (20) 

24^ E^ 2 (21) 
WL = — 2 V^fsm — -- . (22) 



J2\ 7T / ^ v 3 l 2 \ hu c 4 
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In detail, we write the sum in eq. I|18[) as 

oo 

2 y/e-(2n + l)n = (e - tt) 1 / 2 + V(e; 0), (23) 

oo 

ip(s;x) = ( £_ n ~ 2ir\n + x\) 1/2 . 

n— — oo 

Since ip(e; x) is periodic in x, we can use it for the Fourier expansion of eq. I|23(l . and then set 
x = 0. By taking the real part (Re) of eq. (|2"5|) and using eq. (fTTJl . we obtain 

Re{ Eq. ((23I>} = - (' dT(e-T) 1/2 + -Y(-lY f dr(e - t) 1 ' 2 cosz/t, (24) 
^ Jo 77 ~i J° 

where we neglected tt against s. The integral in the first term in eq. (|24ll yields (2/3)e 3 / 2 , 
leading to Wq in eq. I|2Q(I . The integral in the second term can be written after integrating by 
part, and changing the variable [ye — vt = t) as 



1 



cost f U£ sin£ 

sin^e / at — F — cos^e / at — = 
y/t Jo Vi 



We now use asymptotic expansions for ve = x 3> 1: 

' , sint pn cos a; 
dt 

dt 



Vt V 2 
cos t pK sin x 



Vizi) 



Vt V 2 

The second terms in the expansion lead to W-^ in eq. I|21|l . where we used the identity 

_ 12' 

i/=i 

The first terms lead to the oscillatory term W^ sc in eq. (|22H • 

The term Wq, which is independent of B, gives the weight equal to that for a free-electron 
system with no field. The term proportional to B 2 is negative (diamagnetic), and can 
generate a Landau diamagnetic moment. 

The energy E of the 3D system is continuous, and hence the system is manageable or 
soft. In contrast, the energy E of the 2D system is discrete, and hence the system is less 
manageable. This explains the absence of Landau diamagnetism for the 2D system. 

The statistical weight W osc in eq. I|17|l has a sine term. Hence, the density of states, 
J\f = dW/dE, has an oscillatory part of the form 

{2ttE\ n 2 
sin - — , E=- . 

\ nuj c ) 2m* 

If the density of states oscillates violently in the drop of the Fermi distribution function: 
F(E) = \p.Pi E ~^> + l] , the delta-function replacement formula 

-£-*<*-/.). 
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cannot be used. The width of —dFjdE is of the order /c B T. The critical temperature T c below 
which oscillations can be observed is £; B T C ~ T%uu c . For T < T c , we may proceed as follows. 
Let us consider the integral 



/•OO 

1=1 dE F(E) sin (2irE/huj c ) 
Jo 



We introduce a new variable £ = f3(E—fi), and extend the lower limit to — oo (low-temperature 
limit) so that 

1 1 f°° ll/- 00 1 



With the help of the standard integral 

d( 



+ 1 i sinh7ra ' 
we obtain 

7rfc B Tcos(27re F /7kj c ) 

sinh(27r 2 fc B Tm*/?ie J B)' 1 j 

For very low fields the oscillation number in the range fc B T becomes great, and hence the 
sinusoidal contribution must cancel out. This effect is represented by the factor 



[sinh(27r 2 fc B TTO7fr eJ B)] 1 



We calculate the free energy indicated in eq. JHJ using the statistical weight W in eq. I|15(l , 
and obtain 

T = iV M - 2A^e B + 2 A^ fcB T . V til . ^iT^t 5v ( 26 ) 
7T?i 2 tt?i ^-J smh(2ir 2 vk B Tm*/heB) y ' 

where we used the integration formula (|25|l . and took the low-temperature limit except for 
the oscillatory terms. The magnetization A4 can be obtained using eq. |JB}. 
So far, we have not considered the Pauli spin magnetization [7] 

.Mp au ii = 2hIBAT q (e f )/A = 2n^|B/e F , 

with fiB (Bohr magneton) = eh/2m and n (electron number density) = m* e F / 'nTi 2 . Using 
eqs. iJBJ, (J7J and (j^SJl, we obtain the total magnetization M.% t = A^p au ii + Mosc'- 



Mtot = 2n/is 



fc B T m cos (27rep/?ia; c ) 
/i B _B m* sinh (27T 2 /c B Tm* /he B) 



(27) 



The independence of F is contained in the last term in eq. (|26|l . The linear B-dependence 
of the multiplication factor is much stronger than the B-dependence of the alternating series. 
Therefore, the contribution from the £?-derivative of the series is neglected. The variation 
of the statistical weight W is periodic in £? _1 , but it is far from sinusoidal. Only the first 
oscillatory term, v = 1, is important and kept in eq. (|27|l since smh(2ir 2 k B Tm*/heB) > 1. 
The width of dF/dE is finite for a finite T. In this E'-range, many oscillations can occur 
if the field B is made low. We assumed this condition to obtain eq. 127(1 . The magnetic 
susceptibility x is defined by the ratio \ = M./B. 

In conclusion, the 2D system is intrinsically paramagnetic since the Landau's diamagnetic 
term is absent, but the system exhibits a dHvA oscillation. 
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